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Abstract 

In this paper we revisit the smallest-eigenvalue distribution of the Laguerre ensemble by 
presenting in closed form certain integrals obtained previously. With this information we 
compute, using Dyson's continuum approximation, the two-smallest-eigenvalue distribution 
of the Laguerre ensemble. Higher-order contributions to the free energies describing these 
two probabilities in certain scaling limits for (3 ^ 2, which can be interpreted in this context 
as entropic effect, are found. 
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1. Introduction 

In a previous paper [ffl] the smallest-eigenvalue distribution function for the Laguerre Ensemble 
was found via Dyson's || continuum approximation in a certain (double-) scaling limit (to be 
described below). The continuum approximation treats the eigenvalues of an ensemble of N x N 
random matrices by a continuous fluid in the limit large N, with its thermodynamic free energy 
completely characterized as a functional of the fluid density. Although a rigorous proof of the 
general validity of such an approximation is not available, we have confidence in the robustness 
of this approach, since results obtained using it ||, |], ^, ||, |j] compare favourably with those 
obtained from exact analyses, whenever such analyses are available. 

Quite generally, the probability functions which concern us can be computed as the change in 
the free energy when the support of the fluid is perturbed. In most situations the determination 
of the smallest-eigenvalue distribution and other allied quantities can be posed as the following 
problem in equilibrium statistical mechanics: what is probability, prob(-R), of finding a bubble 
of radius R due to thermal fluctuations in the bulk of a fluid at temperature T? This important 
quantity measures the degree of correlation in the fluid. As the creation of a void requires 
energy, such events are relatively rare. Of particular interest is the behaviour of this probability 
when the bubble is large. A physical argument that gives a feel for the answer goes as follows: 
creating a bubble of radius R in d— Euclidean dimensions requires a volume energy oc R d and a 
surface energy oc R d_1 . An application of the Boltzmann principle then suggests that 

prob(ii) ^e~ ClRd - C2Rd ~\ (1) 

where c\ and C2, being respectively the volume and surface energies divided by ksT, are positive 
constants independent of R. Eq.(l) is valid for a bubble with a radius much larger than the 
coherence length. Thus, prob(-R) ~ e - constantx ^ j n Q ne dimension. For a fluid with short-ranged 
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interactions Eq.(l) is expected to hold, but it should fail for a fluid with long-ranged interactions 
as is the case for a charged system, because the volume and surface energies are then no longer 
well defined. That such a consideration is pertinent to random-matrix problems is due to the 
appropriate measure in the space of N x N matrices consistent with symmetries. The invariant 
measure, according to a classification theorem of Dyson Q is 

N 

P(M)dM = exp [-tra(M)] dM = J] \x a - x b f J| exp [-u(x c )] dx c dG, 

l<a<b<N c=l 

where x a are the eigenvalues and dG is the Haar measure of the symmetry group that diagonalizes 
the matrices with orthogonal, unitary and symplectic symmetries, corresponding to (5 = 1,2,4, 
respectively. 

The probability that an interval J (a subset of R) is free of eigenvalues, also known as the 
Gap Formation Probability (GFP), is then 

where J is the compliment of J, JUJ is natural support of the eigenvalues, d/j,(x) := exp[— u(x)]dx 
and 

W(x!, . . . ,x N ) := -(3 lD.\x a -x b \. (3) 

l<a<6<7V 

Insight into this distribution can be gained by considering the x's as the spatial coordinates 
of the particles and u(x) can then be interpreted as a form of "confining" potential holding 
together the logarithmically-repelling charged particles (all carrying charge +1) moving in one- 
dimension. Under the combined influence of the confining and repelling forces the charges will 
attempt to settle down to an optimal configuration with small fluctuations around it. In the 
limit large N we can treat this system as a continuous fluid using techniques from electrostatics 
and thermodynamics. From the above consideration, the GFP can be written as 

E p [J\=e- SF W, (4) 



where 

6F[J] := F[J] — F[J U J], 

is the free energy of N charges confined to region J minus the free energy of N charges in the 
natural support of dfi. As indicated above, the heuristic formula — Eq.(l) — for the GFP fails 
for systems with long-ranged interaction. Instead, an alternative argument peculiar to charged 
systems which provides a qualitatively correct estimate for a large interval can be formulated 
as follows. An appropriate counter charge distribution is introduced externally to neutralize the 
region from which we wish to have the eigenvalues removed. If Fq is the free energy of the 
system with the addition of the counter charge (Q) and Fq = q the free energy of the undisturbed 
system then, according electrostatics, Fq differs from by an amount equal to the charging 

energy, i.e. 

FQ ~ Fq + -^charging, (5) 

where 

Q 2 

-Charging = 7^°<AA 2 [J], C = capacitance 



[the total number of eigenvalues excluded from J] 2 



dx a(x) 
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(6) 



and <t(x) is the eigenvalue/charge density []. Numerous examples supporting Eq.(6) can be found 

in @, ||. A further example that deals with the Laguerre ensemble is given in Appendix A. We 

mention here that a screening theory of the charged fluid gives a physical justification || of 

Eq.(2), although we do not know of a proof of the general validity of this relationship. 

1 Exact Painleve analysis jjj] for large |J| and for /3 = 2, whenever such treatments are available, show that 

there is a sub-leading term ~ lnA^J]. This is not obtainable from the charging energy arguments. However, as 

was shown previously |l| and will be shown later for the Laguerre Ensemble, the logarithmic correction is obtained 

in the continuum approximation. This is due to the fact that the density of the Laguerre Ensemble is very large 

at the "hard edge." 



4 



In this paper we will use Dyson's continuum approximation, so let us first briefly outline 
what that entails. Dyson showed that the functional dependence of the free energy on the charge 
density is given by the functional 

F[a] := f dx u{x)a{x) — — [ dx [ dy <j{x) In \x — y\c(y) + ( 1 — ^\ [ dxa(x)ln[a(x)], (7) 

JK 2 J k J K V 2 / JK 

where the fluid is confined to the region K and the third term on the right-hand side of Eq. (0) 
came out of Dyson's careful analysis and will be referred to as the entropic term. Given this 
functional, then the equilibrium density <r(x) is determined by minimization of F[cr] with respect 
to the (j{x) subject to the normalization constraint 

/ a{x)dx = N, (8) 

JK 

N being the number of charges, which is equal to the rank of the underlying matrices. Using 
Lagrange's method to handle this constrained minimization, which is equivalent to introducing 
a chemical potential, A, shows that the integral equation satisfied by the ground-state density 
is 

u(x) - J dy a(y) In \x - y\ + ( 1 - In a(x) = A. (9) 
Therefore the free energy at equilibrium and with precisely N charges contained in K is 

F[K] = -NA + - f dx u(x) a(x) - ( 1 - ~) [ dx a (x) In a (x), (10) 
2 2 Jk \ 2 / Jk 

the calculation of which for various different cases will form the core of this work. 

This paper is organized as follows: in Section 2 we revisit the Laguerre ensemble and compute 
the probability that the interval [0, a] is free of eigenvalues by utilizing closed form expressions 
of certain integrals that were previously determined approximately and supply in a (double-) 
scaling limit the change in free energy including the entropic (/? ^ 2) contribution; with this 
information we then present in Section 3 the two-smallest-eigenvalue-distribution function for 

5 



the Laguerre ensemble; in section 4 we present our conclusions and discuss the extension of our 
method to the m-eigenvalue problem; details of the calculation of certain key integrals are given 
in the appendicies, along with a note on the charging energy of the Laguerre ensemble. 



2. The Laguerre Potential Revisited 

The confining potential u{x) for the Laguerre Ensemble is 

u(x) = x — alnx, x > and a > —1. (11) 

The change in free energy is computed as 

SF = F[b,a]-F[b,Q], (12) 

where F[b, a] is the free energy where all N charges are supported on a < x < b. This would 
require the determination of the fluid density <r(x) which satisfies the singular integral equation^ 

^=/5P f'dy^yl, a<x<b, (13) 
ax J a x — y 

the general solution of which for the ground state is Q 



ti fJ y ju ti Jay Ju y u y Liy 



For the Laguerre potential the fluid density is therefore 



, v 1 Ib-x 



a a 



(15) 



7T/3 y x — a 

where we require that y/ab > a so as to ensure that <j(x) is positive semi-definite. The normal- 
ization condition now reads 

b — a a ( fa \ 



2/3 p VV b 



This equation is found by differentiating Eq.(9) with respect to x, for (5 = 2. The (3^2 case will be considered 
later. 

3 c/ y/Jb — x)(x — a) (c a constant) solves the homogeneous part of Eq.(|l^), but including this would clearly 
increase the free energy and so it must be excluded from the solution for the ground state density. 



which gives b as a function of N and a via a solution of the following cubic equation (found from 
a slight rearrangement of Eq. ( |i~6| ) ) : 

b 3 - 2{u + a)b 2 + {v + a)b - 4a 2 = 0, := 2/3iV + 2a. (17) 



For sufficiently large N this has a positive discriminant and therefore three real solutions, de- 
terminable by the trigonometric method || . Since it is a simple matter to convince oneself that 
b decreases on increasing a and that the leading order contribution in a should therefore be 



negative it follows that the root which we require is 



2{y + a) 



1 + cos I - cos 



-1 + 



54a 2 a 



2vr 

T 



(18) 



2.1 The Free Energy 

Given b and cx(x) we can now determine the free energy which is given by the sum of the chemical 
potential contribution -F c hem : = 

NA/2 = A and that due to interaction, i.e., 

Fmt ■= \ f <j{x)u{x)dx. (19) 

To compute the chemical potential we first drop the /? ^ 2 correction in Eq.(||) (we will return 
to this term in Section |2.3|) and then evaluate it with x — > b giving 



A = b — a In b 



a (b 
In 



4e 



+ OL 



1 



+ ad - h 



b — a 



(20) 



where I\(x) is tabulated along with the other integrals Ij(x), j = 2 ... 4 in the Table, whilst the 



details of their computations are confined to Appendix B. Using the Table we have 



chem 



V — OL 



b — a\, fb — a\ [a fb — a 

In — - — — aw - In — - — 

2/V4e7 V6V4 



2a tanh 



(21) 



The interaction contribution to the free energy is slightly more complicated; 

„ (b — a) 2 b — a v ^, , s1 a(b — a) , / a 

Fint = TTTn + —T7T i a ~ a ^ b -<*Hb ~ a)] - no h ' 



16/3 4/3 



20 



b — a 



+ o?Hb-a) A ffWJ. , (22) 



2/3 V V b J 2/3 \b-aj ' 
where ibO* 1 ) and ^(x) are in the Table. 

2.2 The Scaling Limit of F[b, a] 
From the expressions for Fchem 

and F int we shall now obtain the free energy in the double scaling 
limit, i.e., the limit where ./V (and therefore v) — > oo and simultaneously a — ► 0, such that the 
combination s := va is finite. Under the above limit the free energy will become a function of 
iV plus a function of s only; other terms of 0(a^ fJ.,S > are discarded. Thus, expanding 
the free energies in a, we get 

F[b,a] ~F (N) + ^ lns + y v^+^, F (iV) ~ -|iV 2 IniV. (23) 

The GFP is thus given by 

exp [-A + 

in agreement with previously obtained results 0]. 

2.3 The Entropic Contribution to the Free Energy for (3^2 

The results we have so far obtained are strictly only true for (3 = 2, as for (3^2 there is an 
extra entropic contribution to the free energy and in the previous sections we have ignored it. 
However, if we assume that the factor (1 — (3/2) is small, then to lowest order in this factor we 
can obtain the total free energy by adding in the term 



ent •" 



1 " 2) f a ^ ln(7{x) dX =( 1 - f ) /«* ( 25 ) 



to F, the density being as per Eq.(|15|), and, with A := 5z^y f and x := 

/ ent :=-ln(vr/3)iV + /« + /(2) ) (26) 



where 



/ <.> := ^( I (itl /iH<fi-^W ,1 - t 



7T/3 Jo 



t + X 



(27) 



and 



7T/3 Jo 



A \ / A 
In 1 



t + x 



(28) 



We start with f^t which can be straightforwardly reduced to 



f (i) 

J cnt 



b — a A(b — a) 



P 



2(3 



[h{x)-h{x)\ 



— a a (a 



(29) 
(30) 



Thus, in the scaling limit, we find that 



J cnt 



v a ( s 1 / 2 

~ T, + In S + O 



P 2/3 



v 



(31) 



(2) / 

Concerning / en t, we first observe that A/x < 1, because of the constraint a < y ab imposed to 



ensure that the density is non-negative, which also makes certain that logarithm in F ent never 

A ( < ?\ 

leads to imaginary terms. Thus, for t > 0, j-^ < 1 and we can therefore evaluate _f en t by 



expanding the logarithm in Eq.(|30D as a power series in In the scaling limit we find that 



n+l 



/3^4« [(n + l)!]' 



1\ 1 
1 + - — , 

nj S 2 



(32) 



the details of the derivation being in Appendix C. We have, upon collecting Eqs.(p3[), ( pl| ) and 



Ep{s) ~ c- 



cxp 



2/3 ^ /3 V* 



q2 . «(2-j?) 

S 2/3 4/3 



1 + 



(33) 



c being a constant. 

To conclude this section we note that from a probability argument we can also obtain the 



smallest-eigenvalue distribution py(s), as it is related to Ep{s) by 



vf{s) 



dE p (s) 
ds 



(34) 



Thus, from Eqs.(|33[) and (|34]), we find that 



P /9 



2a 

2/3 V 7¥ + 



a + 1 



/3 



exp 



s ,2a 2 
20 + JT b 



a 2 a(2-0) 
S 2/3 4/3 



(35) 



3. The Two-Smallest-Eigenvalue Distribution 



(2) 

In this section we will calculate, within the continuum approximation, p^ (x±, X2), the distribu- 
tion function of the two smallest eigenvalues x\ and X2 (x2 > x\), given that the eigenvalues are 
distributed according to the joint probability distribution function 



1 



N 



p( Xl ,...,x N ) = ^l[e- u ^ J] \ x * 



I/ 3 



Z 



(36) 



a=l 



l<b<c<N 
(2), 



where Z is the partition function. To calculate p^ (xi,^) we fix x\ and X2 and integrate 
p(xi, . . . ,xn), with respect to X3, . . . , xm from X2 to 00, then 



pf( Xl ,X2) = N(N „~ 1) \x 1 -X2fe- u ^- u ^x 



N 

n 

a=3 



2Z 

dx a e~ u{ - Xa) \xi - x a \P\x2 - x a 



n \ x * 

3<b<c<N 



(37) 



Shifting the integration variables using x a = t a + X2 and a relabelling the indices, we get 



pf(x 1 ,x 2 )=C\x 1 -X2\ l3 e- u{xi) - u( - X2) x 



N-2 r roo 

IT / dt a e~ u ^ +x ^ \ Xl -X2-t a \P\t a \P 



n h-tcf, 

Kb<c<N-2 



(38) 



where C = — ^jz~~ ~- The multiple integral after the trivial factors will be evaluated using the 
continuum approximation described above. This is particularly interesting since the logarithmic 
(one-particle) factors already appear here. We expect, from the results obtained in Section 2 
that there will be an analogous logarithmic correction to the free energy. It is convenient to next 
transform to the centre-of-mass and difference coordinates, e = X2 — x\{> 0) and R = {X1+X2) /2 
respectively, and then take pp(e, R) = p*\xi, x^)- Then 

pf(e, R) = Ce 13 exp -u (r - |) - u \R + | 
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'N-2 



n dt ° 



a=l 



e . u{ta+R+e/2){ta{ta + £))P jj % _ t ^ 

Kb<c<N-2 



(39) 



The multiple integral in Eq.(|39|) can be interpreted as the partition function for a charged fluid 
of N — 2 of particles in the effective external potential u e g(t; e, R) defined by 



u eS (t;e,R) =u{t + R + e/2) - 0hx[t(t + e)}. 



(40) 



For reasons which we will explain later, the potential with a = is much simpler than the 

(2) 

case, so we will first calculate the asymptotic expansion of for a = 0. We therefore have 



pf(e,R) =Ce /3 e" 

(N-2 

i n dx ^ 

l a=l 



2R S 



3 -(a; a +R+f ) 



n \ xb ~ Xc 

Kb<c<N-2 



The analogous integral equation for electrostatic equilibrium is, 



A = x + R -\ /3 In x — ln(x + e) — / dya{y) In |x — y \ 

2 Jo 



where A is the chemical potential, or 



i-/>r± + ; 



a; x + e 



/3P Tdy 
Jo 



<K2/) 



following the method outlined in Section 2. The solution of Eq.(S3) is 



<t{x) 



1 lb — x 



tt0 y x 



/9 



x + e V ^ + e 



, < x < b. 



With the normalization Jq dxa(x) = N — 2, we find 



N -1 = — + 



20 V b + e 



which is equivalent to the cubic equation 



X 3 - 2(y + e)X 2 + (u + efx - 4/3 2 e = 0, X:=6 + e, v:=20(N-l). 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



4 Recognising that we could obtain an exact solution for b prompted us to look for closed-form expressions for 
the integrals and hence to re-examine the level-spacing distribution, the results for which we presented earlier. 
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Again the discriminant of Eq.(46) is positive for sufficiently large N and so we have three real 
roots. The appropriate solution is, 



2{v + e) 



1 + cos ( - cos" 1 (-1 + 54/? 2 e/(^ + ^) 3 ) - — 
3 ^ ^3 



e. 



(47) 



The analogous free energy is 



where 



~ 4/3 l ' 



f f b 

Fint-=-x dxa(x)u efi (x;E,R), 
2 Jo 



(48) 



(49) 



is the appropriate interaction energy. 



3.1 The Chemical Potential A 



To compute the chemical potential we take x = b in Eq.(|4"2]), then 



A= l ^b + R+ £ --PHb(b + e)]- b -^ + pinb 



1 



l + x 



l + x 



h(x), (50) 



defining x := | . In the scaling limits v —* oo and e, R —* we find two suitable scaled variables: 
Si := ve and 52 := + e/2) (for reasons which will be made clear below), yielding 



-Fchem ~ -— In ( — \ - - ln(4e) + -lnu 
8/3 \Ae) 4 V ' 2 



2 4/3 I v 



1/2' 



(51) 



as the asymptotic expansion for F c h em . 



3.2 The Interaction Energy 



In order to simplify the notation we first split up Fi nt into the sum of two terms, Fi and F\^^ 
where 

1 fi fT~f / / g /3 \ r e i 

1 " i r> ' 1 " ' (52) 



2k (3 Jo 



f 



b + et + x 



bt + R + --2dlnb 
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and 

F ^-k!l dt ^( h -^tT-^ lnlt(t+x)] - (53> 

Evaluating these integrals and expanding the results we find that: 



e 



2 " ln " ^ +§ + <>(-) (54) 



16/3 2 2 4/3 V f 



and 



(55) 



F ln = M2)6-^ 3 (x)-/31n6(y^-l) + ^/ 2 (x) + ^y^/ 4 (x) 
~ ^ (Vln2- iln(4e)^ i/- 2/3 In i/ + /3 In Si - VSi + O ■ 
The asymptotic expansion of F{ nt is thus 

F int ~ J F; nt (i,) + /31n5 1 -^y5 1 "+^. (56) 



From Eqs.fliq), (|5l|) and d5q) we therefore find that asymptotically 



F(u;S l ,S 2 )^F°(u)+pin(S l )-2^S' 1 + ^. (57) 



3.3 The Entropic Contribution to the Free Energy 



As explained in section |2.3| when (5 ^ 2 there are additional terms in the free energy, the lowest 



order correction resembling an entropic contribution and is given by Eq. (|25j) . We can once again 
split up F cnt as per Eq. (]26|) , with f^ t and f^ t as before, except for that the prefactor is changed 



according to £-# — ► \ and we now take A = £ A /tt - - It follows that we now have 



and 



/«~^-lni/ + iln5i + ^| (5N) 



^[(n + W ^ 
13 



the last line of which is obtained in Appendix C. 
Thus, for a = 0, but (3 ^ 0, the free energy is 



F(v;S u S 2 ) ~ F M + 
1 » 



lnSi-2VSi + -±- + 



2/3 



E 



(2n)! 







3.4 The Distribution for a ^ 

With the effective potential given by Eq.(f40|) and the charge density is given by 



(60) 



1 b-x 



<T(X) 



1 



a 



R + l 



x + R+ | y b + R + § i + eVHe 



/37r y a; 

from which it follows that the normalization condition now reads 

b 



(61) 



N — 1 



+ 



a 



1 



R+ £ 2 



b+R+^r 



(62) 



2/3 \jb + e (3 

where, for ease of notation we have introduced r := R + 1. In contrast to the above calculations, 



we cannot re-express Eq.(62) in terms of a cubic equation, but only in terms of an 8t/i-order 
equation, which is why we previously described this problem as being more difficult than the 



others. However, if we define v := 2[3{N — 1 + then we can rewrite Eq.(32) as 
which we can solve for 6 as a power series in u, the result of the calculation being 



(63) 



b = v-2 



an + (3e2 «2 + /?£2 — 2(a + /3)(ara + /3e2 ) 



i 

Z/2 



+ 



3 
U2 



+ 



o(V§) 



(64) 



which for future reference we write in the compact form b = v + 5b. 

The chemical potential and interaction energy are calculated in precisely the same manner 
as before, so we evaluate the former at x = b and split up the latter into the sum of F± and F\ n . 
Introducing the further definitions 



n 



b + r 



and £ 



b + e 



(65) 
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then the chemical potential is 

A = ^]n(4e)b + r-a]n(b + r)-p\n[b(b + e)]+]nb(a(l-K)+/3{l-£)-^ + 



i ln(4e)6 + r - (/3 + ^ + n K + .j£)]ub + In ( ■ 1 ) [n 7? + . tf] 



2(atanh~ 1 :K + /3tanlr 1 £). 
Asymptotically we therefore have 

F c hem ~ ^chcm( l/ ) 



Si Oi\JT>2 S 2 



2 2/3 4/3 

As for the interaction energy, it is simple to arrive at 

b 2 b . m „„, aril - 11) + Ml - E) N-2. , 
Fi = — - —(atR, + p£) ; / v - + (r - (a + 2/3) In 6), 



16/3 4/3 



2/3 



and thus 



/ 5i a-v/52 52 
2 2/T + 4/3' 



In full, the expression for F\ n is 



Fm = - (l + |) ln(6) Q - a(l - fc) - 0(1 - £)) + 



2/3 \2 



a 
2 



Kh 



r 
b 
e r 
Vb 



s 








)) 


+ 






'( 


r 






b 


I) 



cr (r 



+ 



which yields the asymptotic expansion 



(66) 



(67) 



(68) 



(70) 



(71) 



F ln ~ F° n (v) -y/S~i + (31n(Si) - ^ + \ (l + ~) ln(S 2 ) + 2a ln(v^i" + v^)- (72) 
So, from Eqs.(p8|), ( |70| ) and (f72[) , the asymptotic expansion of the free energy for non-zero 



a is 



F(u; Si,S 2 ) ~ F u (^)+/31n(5i)-2 K/ft + Wft + 7775 + 



o 



(3 V V 2/3 
1 + ^ ln(5 2 ) + 2a ln( v 7 ^ + v 7 ^)- 



(73) 
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3.5 The Distribution of the Two-Smallest Eigenvalues 



To finally obtain the asymptotic result for the two-smallest-eigenvalue distribution we first intro- 
duce Si := vxi, i.e., the double-scaled variables of the two smallest eigenvalues. Next, we define 

(2) (2) 

Pq (si, S2) = lim i ,_ >0O pg (si/u, sijv), which is the double-scaled asymptotic of the distribution. 
From Eq.(|39|) and the observation that S\ = s 2 — s\ and S2 = s 2 , we therefore have 



pf\a; Sl ,s 2 ) =Ce^^ S2 - Sl '>- F( - u ' S2 - sl ' S2 \ 



Explicitly, the asymptotic distribution with a = 0, but (3 7^ 2, is 



(74) 



Pg 2 \0;s 1 ,s 2 ) = C u exp ~ (l - ln(s 2 - si) + 2y/s 2 - si 

l\ ^ (2n)\ ( 



S2_ 



E 



13 



where C v represents all the Sj-independent terms. 

Similarly, with and (3 = 2, from Eqs.(73) and (74), we then have 

P 2 (2) (a;si,s 2 ) = C^exp + 2^s 2 - si + a (y/s~2~ - i (l + ~ \ ln(s 2 )^ - 
2a In (j~s~2 - si + v^2^ 

3.6 A Comparison of the Continuum Approach With Exact Results 



(75) 



(76) 



Exact results for the distribution of the two smallest eigenvalues have been obtained by Forrester 
and Huges (FH) [O] for the unitary Laguerre ensemble with a = t, where I is either zero or a 
positive integer. We are thus able to check our results and will do so for the cases a = and 



a = 1. In order to proceed we first note that in terms of the notation in [12|, we have a = 1 and 
m = N. Secondly, the variables Sj introduced in the previous sub-section are, in the limit of 
large N, the same as the Sj used in [12|; we are now ready to check our results. FH have shown 
that an exact expression for the scaled asymptotic distribution is 

P 2 {2) {l; Sl ,s 2 ) = 2- 4 e - f r +£ln ( S2 / Sl )x 
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det 



j = !,...,£ 
k = 1, ...,£ + 2 



(77) 



3 = 1,2 
k = l,...,e + 2 

which requires us to compute the determinant of an (£ + 2) x (I + 2) matrix of modified Bessel 
functions of the first kind, I v (x). Note that these are not to be confused with those in the Table. 
For a = and a = 1, we thus find that 



P (2) (0;ai,a 2 ) 



_£2 

e 4 
16 



h {\fs~2 ~ Slf - h (VS2 - Si) h {\fs 2 - Si) 



(78) 



and 



P (2 \l;si,s 2 ) = 



e 4 



+ln(s2/si) 
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— — ) h{s 2 ) \h {\fs~2 ~ Si) 2 ~ Iq (y/s 2 ~ Si) I 2 {\J~S~2 ~ Si) 

S2 



I S2 ~ Si 
S2 



h{y/S2) [h WS2 - Si) h (VS2 - «l) - ^0 {\fs~2 ~ Si) I 3 (V«2 = Si)] + 



h{\fs~2) h {\fs~2 - si) 2 - h (\fs 2 - Si) I 3 (VS2 - Si)] } 



(79) 



respectively. We claim that the Coulomb-fluid gives the asymptotic limit of the given distribu- 
tion, i.e., the distribution in the limit of y/s 2 — si, y/s2 3> 1. In the limit of large x, the modified 
Bessel functions have the asymptotic expansion 



L(X) 



V2 



1 



TTX 



8x 



+ 0(x- 2 ) 



(80) 



The asymptotic expansion of the given distribution is therefore 



P 2 (2) (0;si, S2 ) 



; -s 2 /4+2Vs2-si 

32^ 



S2 



^ + 0((S2-S 1 )- 2 ) 



(81) 
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for a = 0, whilst for a = 1 we have 



P 2 (2) (l; 8l) s 2 ) ~ ^4e 2 v / ^+v^ s -f [i + . . .] . ( 8 2) 
16(2vr)2 

^From the Coulomb-fluid approximation we get 

P 2 (2) (0; si , s 2 ) ~ e - f +2v^T ( 83 ) 

and 

P 2 (2) (l;Sl,S 2 ) ~ e"f + 2 V^2- S i+v^ x e -|ln( S2 )-21n(V S2 -s 1+v ^)_ ^ 

Comparison of Eqs.(|8~I|) and (83) and Eqs.(^) and ( |8"4| ) reveals that the continuum approxima- 
tion successfully captures the dominant term in the asymptotic expansion of the Bessel functions, 
i.e., the term e~ x , although the logarithmic factors are not precisely reproduced. 



4. Conclusions 

In this paper we have used the continuum approximation to calculate the GFP, and hence the 
smallest-eigenvalue distribution, and the two-smallest-eigenvalue distribution of the Laguerre 
ensemble of random matrices. 

The GFP has been the subject of a previous paper Q, but here we presented a calculation in 
which key intermediate results were given in closed form, in contrast to the simple power series 
which we had formerly obtained. 

The two-smallest-eigenvalue distribution was new and demonstrates how it is possible to 
extend the continuum approximation to deal with higher-order correlation functions. The two 
essential differences in the calculation being: (i) the normalisation condition is that there should 
be N — 2 charges in the interval and (ii) the set of the smallest eigenvalues modify the potential 
in which the remaining charges are placed, thereby leading to an effective potential instead of the 
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plain Laguerre potential. Exact results for this distribution are available for (3 = 2 and a = or 
a positive integer and we have shown that the continuum approximation gives the leading order 
terms in the asymptotic expansion of this distribution for a = and 1. We therefore feel justified 
in claiming that our result is not only correctly gives the dominant terms for a = 2, 3, . . ., but 
that the leading-order factor in the asymptotic expansion of this distribution for all a and (3 is 
given by 

(85) 



Pi 2 \a;s 1 ,s 2 ) ~ exp 



■S2 



— + 2y/s 2 - Si + Oty/S^ 





Thus, although the Coulomb-fluid is limited to the dominant factor, the method is somewhat 



simpler than the known method of deriving the distribution exactly [12] and also has the ad- 
vantage of being applicable, we believe, to the full range of parameters of interest. 

Since we have shown that the Coulomb-fluid approximation can be applied rather successfully 
to problems regarding the smallest and the two-smallest eigenvalues, it is only natural to ask 
whether we can use it examine the m-smallest-eigenvalue distribution. If x\ < x% < ■ • • < x m 
denote the m smallest eigenvalues, then, from Eq.p^), the m-smallest-eigenvalue distribution 
can be written in the form 

p^( Xl ,...,x m )=C^ J] \x a -x b fe-^<^vf\ (86) 

l<a<b<m 

where is the normalization constant and 



X>( m ) : = ( dx m+ i ... I dx N 



Y[ \x b -x c f Y[ \xd-x 

l<b<m,m+l<c<N m+l<d<e<N 



(87) 



Introducing the change of variable t a - m = x a — x m , we have thus transformed the problem of 
calculating p^ 1 to the problem of determining , where 



Kb<c<N-m 
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the effective potential now being 

m 

u cS (t;xi, . . .,x m ) :=u(t + x m ) - (3^\n\t + x m - x a \ . (89) 

a=l 

We propose that X>l m ^ can be calculated asymptotically using the Coulomb-fluid approximation 
with the normalization being set equal to N — m; we will not persue this problem any further 
here. 
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Appendix A. The Charging Energy of the Laguerre Ensemble 

The integral equation for the density, Eq. (9), for (5 = 2 and with the specialization to the 

Laguerre Ensemble (u(x) = x — a In x) reads, 

,6 

x — alnx — (3 1 dy a(y) In \x — y\ = A, < x < b. (Al) 
Jo 



Eq.(Al) in addition to regular solutions also admits solutions that describe point charges. In- 
cluding these would produce infinite self-energy thus causing the total free energy to diverge. 
Bearing in mind that the GFP is expressed as the difference of two free energies, see Eq. (4), the 
self-energy divergences cancel. With the point charge at the origin superimposed on the regular 
solution, Eq. (Al) has the following solution: 



a „ , . 1 b — x , , . , 

a(x) = --8+(x) + — \ , 0<x<b. (A2) 

p 7rp y x 

Observe that the point charge carries weight of order unity while the normalization condition 
Jq a(x)dx = N, is almost exhausted by the regular distribution. We find from the normalization 
condition, b = 2(3N + a. A simple calculation shows that the number of eigenvalues excluded 
from the interval, [0, a], (a < b) is 



f a a 1 

Af[a] = J q a(x)dx = - Tp + - p 



a(b — a) + 6 tan 



-i 



(A3) 



b — a 

In the scaling limit b — > oo and a — > such that ab = s is finite, we find 

^i~-i + ^ +0 (S)> < A4 > 

and 

Ep{s) ~ exp [— ois + a,2\/s ] , (A5) 

where a\ and 0,2 are constants independent of s. Eq. (A5) is in qualitative agreement with Eq. 
(20). As noted in the introduction the charging energy argument does not supply the logarithmic 
correction. 
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Appendix B. Evaluation of the Integrals 



The integrals Ii, I2 and Z4 are all evaluated by the same method, whilst I3 is evaluated by ele- 
mentary means and is given later. Firstly, we use the identity /i A In h = d\ (where d\ = ^ 
which for the integrands of interest yields an integral proportional to the Gauss hypergeometric 
function |fTo(|. More precisely, we have: 



h(x) 
h(x) 



7T 



1 



2 V x + 1 



x X 2*1 f-A,-;2; 

2 x 



A=-l 



(Bl) 
(B2) 



and 



J 4 (a;) = — <9 A 
7rx 



(B3) 



Next, we apply the transformation formulas fTT|| which send 2-Fi(-> -j •> to 2-Fi(., ., ., x) to get 
a convergent hypergeometric series. Since this step and those which follow do not essentially 
differ for the three integrals, we only display the calculation for l2{x). The transformation 



formula used is Eq.(9.1322) of [10], from which 



h.(x] 



1 + x 



+ ln(4x) 



(B4) 



where 



and 



h{x) := d x 2 Fi (l - A, - A; | - A; -x) 



A=0 



i 2 (x) := dx 



:l !-x) A+ ^ 2 Fi (a,A + 1;A + -;-x 



A=0 



(B4) 



(B6) 



To compute i\ and 12 we use the series representation of 2-^1 1 which is differentiated to give 



k(x) 



(B7) 
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and 



ln(l + x) + ^B -,n )(-xr 



n=l 



(B8) 



i 2 (x) = y/1 + X 

Using the integral representation, 

B(x,y)= f 1 dtf-^l-ty- 1 , 5ftx, Ky>0, 
Jo 

of the Beta function, followed by exchanging the order of integration and summation, we find 
with the aid of En=o(- xt ) n = (1 + xt)' 1 , Ml < 1, 

h(x) = x f'dt^j, (B9) 

JO 1 + XT 

and 



i 2 (x) = a/1 + 



x 



ln(l + x) — x 

J o 



dt 



y/l-t(l + xt) 

The integrals in Eqs. (B9) and (BIO) are expressed in terms of elementary functions, 



(BIO) 



h(x) = -2 + 2 



X tanh 1 



1 + x' 



and 



i 2 (x) = vT+~x 
and thus the result for -^(a?) follows. 



ln(l + x) - 2 



1 + x 



tanh 1 



1 + x 



(Bll) 



(B12) 



To determine I%(x), we first compute I' 3 (x). The t— integration is then elementary giving 



r 3 (x) = -i + 



l + x 



(B13) 



We find by integrating Eq. (B13) and supplying a constant, 



h{x) = J 3 (0) + £ \j^f- d y = "^y^ ~ x + + + In (v^ + v^T^) • (B14) 

The function I§(x,y) requires a little more effort to obtain its closed form. We begin by 
writing the factor ln(£ + x) = ln(t) + P J X duj (u + t) and thereby obtain 

h(x,y) = h{y) + P / -J— —[——-— — 

Jo u-yJo ttV t \t + y t + u 



(B5) 
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The ^-integration is easily performed and leads to 



h(x,y) = h(y) + P f 
Jo 



du I 1 + y 1 + u 



u-y \ V y 



(B6) 



Given that 



7 

Jo 



p ^/T+, =lii 
u — y V u 



1 + 2x + 2 



1 + x 



+ 



ln(x - y) + In ( -x - y - 2xy - 2^x(l + x)^Jy{l + y) 



(B7) 



then the result for 1$ shown in the Table follows easily. 
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Appendix C. Entropic Contribution to the Free Energy 



(2) 

The entropic term for the smallest eigenvalue distribution is given by 



f (2) 
J ent 



b - a f 1 dt fT^t 



(3 Jo vr 



t 



1 



A 



t + x 



In 1 



A 



t + x)' 



(CI) 



where A := -^y^. We could evaluate this in closed form, but this would require a method 
different to that in Appendix B and as we are most interested in the scaling-limit form it is more 
sensible to focus our attention on that. To this end we proceed by first expanding the logarithm 



in Eq. (CI) as a Taylor series in since \fah > a (the positivity condition on the density) 



implies that < 1 for t > 0. Doing so gives, 



f (2) = I 
J ent p 



-91 + J2 



9n+l 

^ n(n + 1) 



. 9i = a[ x j--l) , 



(C2) 



where 



dt I -t ( 1 



V Jo f V t \t + X 



= d\ 



/Zf(-A) 



n— 1 



d 



b (n — 1)! \dx 



n-l 



1 + X 



, n > 2. (C3) 



Writing ((i/tix)" 1 ^! + % = «m(l + x)a m and (d/dx) m x 2 = c? m x 2 m , then an application of 
the Liebnitz's rule gives 



™ n-l , 

^ = (-lrv (?) E ,r bhp aP ~ n ~ l 



(C4) 



In the scaling limit we see that the only significant term in Eq.(C4) is the p = one. Since 
u = 1 (trivially) and d m = (-l) m (2m)!/(2 2m m!), then 



9n 



a n (2n-2)! 1 

2 2n-2[ (n _ 1)! ]2 (o6)2f 



r + O 



(C5) 



and therefore 



/2~-a + ^E 



a n+1 {2n)\ 



/3 ^ 4™ [(n + !)!]< 



1\ 1 
1 + - — . 

n/ S 2 



(C6) 
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The analogous contribution in the two-smallest eigenvalue problem only differs from the above 
in that the prefactor is ^, as opposed to x = e/b and A = f y^^p- We therefore now have, 

,».._, i f_W (AY 
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Function 


Integral Representation 


Closed Form 


h(x) 


fl dt ll-t ln(l-t) 

JO TT V t t+X 


2 In 2 + ^±±5 [ln(l + x) - 2 tanh" 1 


h(x) 


1 x rl dt jl-t ln(t+x) 

V 1+X JO 7T V t t+X 


2 ^ 2^ + ln[4x(l + x)] - 2 [l + ^ tanh" 1 ^ 


hix) 




_ x + ^(1 + x ) + l n ^ x + VI + x] 


h(x) 


rl dt ll-t lni 

JO TT V t t+X 


2 In 2 + ^±±£ [in x - 2 tanh" 1 ^fj] 


h(x,y) 


fl dt ll-t \n(t+x) 

JO TT V « (t+jO 


2 l n 2 _ 2v /I±^ tanh -i atanh- 1 ^f x + 
^^\n[x + y + 2xy + 2y/x{\ + x) y/y(l + y)] 



Table 1: Table of Integrals. Note that all the resuls are for x > 0. 
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